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ABSTRACT 


11-e  claencal  theory  of  the  scattering  of  a  plane  electromagnetic  wave  by  a  sphere  is  re¬ 
viewed,  and  the  difficulties  in  getting  numerical  answers  from  this  formal  solution  are  dis¬ 
cussed.  It  is  shown  how  the  computations  may  be  simplified  by  using  suitably  defined 
logarithmic  derivative  functions  and  the  recurrence  formulas  due  to  Infeld.  By  this  method, 
numerical  answers  for  the  scattering  amplitude  coefficients  of  any  order  can  be  computed 
exactly,  even  if  the  index  of  refraction  is  complex. 

The  technique  mentioned  above  has  been  used  to  determine  the  scattering  amplitude 
coefficients  and  the  back -scattering  cross  section  for  one  special  case  involving  water  spheres 
with  siaeft  comparable  to  the  wavelength.  The  theoretical  results  are  compared  with  those 
obtained  experimentally. 

A  rigorous  solution  is  also  given  for  the  scattering  of  a  plane  electromagnetic  wave  from 
two  concentric  spheres  of  different  dielectric  constant.  This  problem  is  formulated  in  a 
manner  similar  to  that  for  a  single  sphere,  and  the  scattering  amplitude  coefficients  are  ex¬ 
pressed  in  terms  of  spherical  Bessel  functions  and  the  logarithmic  derivative  functions. 
The  application  to  a  particular  physical  problem  is  indicated. 
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Symbol 

Same 

( 'nits 

a 

Radius  of  the  sphere 

meter 

b.’ 

Scattering  amplitude  coefficient* 

dimensionless 

b 

Radius  of  the  spherical  shell 

meter 

B 

Fundamental  magnetic  vector 

volt-second  square  meter 

D.(x) 

General  Bessel  density  function 

E 

Fundamental  electric  vector 

volt  meter 

//.*(*) 

llanltel  function  of  the  second  kind 

Roman  letters  indicate  unit  vector*  along  the  x,  v  and  x  axes 

/ 

Square  root  of  minus  one 

J.(x) 

Bessel  function  of  the  first  Lind 

k 

Wave  number  of  space 

radian  meter 

b, 

Complex  wave  number  of  the  sphere 

radian  meter 

k» 

Complex  wave  number  of  the  spherical  shell 

radian /meter 

1,  m,  n 

Orthogonal  spherical  vector  wave  functions 

m,  n 

Integer  indices 

J V 

Complex  index  of  refraction  of  the  sphere 

dimensionless 

Nt 

Complex  index  of  refraction  of  the  sphere  (Section  •>) 

dimensionless 

Nt 

Complex  index  of  refraction  of  the  spherical  shell 

dimensionlnw 

Associated  Legendre  polynomial 

P. 

Total  absorbed  power 

wait 

P. 

Total  scattered  power 

wall 

P, 

Total  power  removed  from  the  incident  wave 

watt 

ft 

Total  absorption  cross  section 

square  meter 

ft 

Total  scattering  cross  section 

square  meter 

ft 

Total  attenuation  cross  section 

square  meter 

r 

Spherical  coordinate 

meter 

*.(*),  S,(*) 

Besael  density  functions 

s 

Complex  Poynting  vector 

watt/square  meter 

f. 

Characteristic  velocity  of  space 

meier/setxmd 

d 

Complex  characteristic  velocity  of  the  sphere 

meter/second 

f» 

Complex  characteristic  velocity  of  the  spherical  sliell 

meter /second 

w< 

Power  density  incident  upon  the  scatterer 

watt/ square  meter 

w. 

Power  density  of  the  far  sane  scattered  field  in  direction  of  soum 

x,y,t 

Cartesian  coordinates 

meter 

*.*(*) 

General  spherical  Bessel  function 

Z.(x) 

General  cylindrical  Bessel  function 

a 

Electrical  siae  of  the  sphere 

radian 

p 

Electrical  siae  of  the  spherical  shell 

radian 

<.(*) 

General  logarithmic  derivative  function 

«• 

Dielectric  constant  of  space 

farad  meter 
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C _ 1  -1 

jywmn 

Name 

Units 

«# 

Rest  effective  absolute  dielectric  coostant 

farad/meter 

r. 

Characteristic  impedance  of  space 

ohm 

Spherical  coordinates 

radian 

X 

Wavelength 

meter 

M* 

Permeability  of  space 

henry /me  ter 

r» 

Magnetic  constant  of  space 

meter  /henry 

f 

Complex  dielectric  factor 

farad/meter 

li 

Complex  dielectric  factor  of  the  sphere  (Section  5) 

farad /me  ter 

ft 

Complex  dielectric  factor  of  the  spherical  shell 

farad/ meter 

tr 

Relative  complex  dielectric  factor 

dimensionless 

p.W.  »«(*) 

Logarithmic  derivative  functions 

*. 

Real  effective  conductivity 

mho 

9 

Back-scattering  cross  section 

square  meter 

0 

Solid  angle 

steradian 

u 

Angular  frequency 

radian/second 

Subscript  i 

Refers  to  the  incident  field 

Subscript  s 

Refers  to  the  secondary  field  outside  the  scatterer 

Subscript  t 

Refers  to  the  secondary  field  inside  the  scatterer 

e* 

Asaumed  time  dependence 
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BACK-SCATTERING  OF  ELECTROMAGNETIC  WAVES 
FROM  SPHERES  AND  SPHERICAL  SHELLS 

1.  INTRODUCTION 

When  any  object  ii  placed  in  the  field  of  an  electromagnetic  wave  it  removes  energy  from  the  field. 
In  general,  some  of  the  energy  is  dissipated  intmiailv  as  heat,  and  some  of  the  energ\  is  reradiatrd  to  pro, 
duce  a  secondary  or  scattered  field.  One  of  the  important  problems  of  electrodynamics  is  to  determine 
the  amount  of  energy  absorbed  and  the  amount  and  distribution  of  the  scattered  energy  for  a  particular 
scattering  object  illuminated  under  specified  conditions.  One  of  the  simplest  problems  of  this  tyjvc  to 
visualise  is  that  of  a  plane  wave  incident  upon  a  spherical  particle.  Although  the  geometry  of  this  problem 
is  simple,  the  problem  itself  is  very  important  and  has  received  considerable  attention  from  investigators 
throughout  the  years. 

Pioneering  in  this  field,  Tyndall  (1869a)  showed  experimentally  the  relation  of  the  rise  of  suspended 
spherical  particles  to  the  intensity  and  polarisation  of  scattered  light.  Although  he  Has  mainly  interested 
in  studying  the  chemical  action  of  light,  he  discussed  the  problem  of  the  color  and  polarisation  of  sratlemi 
light  in  several  additional  papers  (1869b,  1870).  Shortly  after  Tyndall's  work.  Lord  Rayleigh  (1871,  1881, 
1897,  1899)  published  a  series  of  papers  which  became  the  basis  for  lie-  classical  theory  of  scattering,  lie 
handled  the  case  of  particles  very  much  smaller  than  the  wavelength  of  light  and  established  his  famous  re, 
suh  that  under  these  conditions  the  scattering  is  inversely  proportional  to  the  fourth  |s>wer  of  tlir  wavelength. 

The  first  thorough  treatment  for  a  spherical  particle  of  any  siae  and  any  electrical  properties  was  givrn 
by  Mie  (1908).  Ilis  solution  is  an  infinite  series  of  spherical-mode  functions  with  amplitude  eurflieirnls 
that  are  determined  from  the  boundary  conditions  at  the  surface  of  the  *|dvrrr.  Kadi  of  these  s|vlieriral- 
mode  functions  itself  involves  infinite  series,  so  that  in  the  general  case  of  a  large  s|dierr  with  a  com|de\ 
index  of  refraction  the  problem  of  getting  numerical  answers  is  quite  involved.  However,  the  Mie  |>a|>rr 
still  remains  the  fundamental  contribution  to  the  subject. 

Since  publication  of  the  paper  by  Mie.  additional  investigation*  have  lieen  msdc  by  Delvyr  (1909), 
Gans  and  llappel  (1909),  Bromwich  (1920),  Ray  (1921),  Jobel  (1925).  Illomer  (1925.  1926).  Stratton 
et  al.  (1930,  1931,  1941),  Trinks  (1935),  Engelhard  and  Fries*  (1937).  Ryde  et  al.  (1911.  1911,  1915,  1916), 
Ruedy  (1941,  1943,  1944),  LaMer  e?  al.  (1943,  1946.  1948).  Brilhmin  (1913.  1911,  1919).  L.  Goldstein 
(1945),  Sinclair  (1947),  Houghton  el  al.  (1949),  Klotxbaugh  and  Duckett  (1919)  and  Olivers.  The  com¬ 
plete  theory  is  given  concisely  by  Stratton  (1941)  and  L.  Goldstein  (1945).  ’flits  will  be  reviewed  in  Sec¬ 
tion  3. 

It  is  not  the  present  purpose  to  review  all  the  papers  dealing  with  live  scattering  problem.  Instead, 
the  discussion  will  be  restricted  to  one  particular  aspect  of  the  (imltlrm. 

With  the  development  of  microwave  radar  during  World  War  II,  it  was  found  that,  five  sufficiently 
small  wavelengths,  rain  could  produce  an  appreciable  echo  ami  substantial  attenuation.  The  echo  plvr- 
nomenon,  was  important  for  several  reasons.  For  the  usual  lartiral  use  of  the  rqui|vment.  it  was  important 
to  distinguish  between  atmospheric  reflections  and  operational  targets.  Mnrrovrr.  as  a  meteorological 
tool,  the  observance  of  echoes  was  very  useful  in  mapping  the  rain  areas  and  in  allowing  their  movements. 
The  problem  of  attenuation  was  important  to  both  radar  and  rommuniralion.  It  is  not  surprising,  liven, 
that  these  problems  received  considerable  attention. 
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To  investigate  theoretically  the  microwave  reflection  from  rain,  it  is  necessary  to  assume  that  the 
urops  are  randomly  distributed  in  space  and  that  the  mutual  interaction  between  drops  is  negligible.  The 
first  assumption  appears  to  be  valid  without  further  comment.  'I  lie  second  assumption  is  based  on  Trink’s 
(1935)  analysis  that  for  Rayleigh  scattering  ihe  mutual  interaction  is  negligible  for  sphere  spacings  greater 
than  two  or  three  sphere  diameters,  plus  the  fact  that  in  actual  rain  the  average  spacing  between  drops  is 
many  times  this  value  (Palmer  (1949)).  Under  these  conditions  the  problem  can  be  divided  into  two 
parts:  (1)  finding  the  drop-sire  distribution  and  (2)  finding  the  reflection  from  a  single  drop.  The  drop- 
size  distribution  is  strictly  a  meteorological  problem;  the  reflection  from  a  single  drop  ia  a  problem  in  electro¬ 
dynamics. 

The  first  thorough  investigation  of  the  reflection  and  attenuation  effects  from  rain  was  made  by  Ryde 
et  al.  (1941,  1944,  1945,  1946).  llis  method  was  to  compute  the  numerical  results  for  a  single  drop  using 
the  Mie  theory  and  to  use  the  assumptions  listed  above  to  make  application  to  the  case  of  many  drops. 
The  main  difficulty  lav  in  getting  the  numerical  results  from  the  Mie  theory.  Despite  the  tremendous 
labor  involved,  however,  he  computed  numerical  results  for  a  variety  of  cases.  For  radars  of  longer  wave¬ 
length,  it  was  found  that  the  Rayleigh  theory  was  adequate,  but  for  radars  of  shorter  wavelength,  a  more 
evact  theory  was  needed.  The  work  of  Ry  de  is  discussed  briefly  by  L.  Goldstein  (1945),  who  has  extended 
it  somewhat,  using  more  detailed  meteorological  data.  Recently  a  table  of  scattering  functions  hat  been 
prepared  by  Lowan  (1949)  and  these  results  have  been  used  by  Haddock  (1947)  in  a  further  contribution  to 
this  probier 

In  all  cases,  a  major  difficulty  has  been  the  evaluation  of  the  reflection  from  a  single  sphere.  It  was 
an  interest  in  this  problem  that  prompted  part  of  the  present  investigation. 

In  Section  2,  the  common  absorption  and  scattering  parameters  are  defined.  Particular  attention  ia 
given  to  the  back -scattering  cross  section,  which  is  the  main  subject  of  this  report. 

The  classical  theory  of  scattering  from  a  sphere  is  reviewed  in  Section  3.  The  difficulties  in  getting 
numerical  answers  are  discussed,  and  it  is  shown  how  the  computations  may  be  simplified  by  using  log¬ 
arithmic  derivative  functions  and  the  recurrence  formulas  due  to  Infeld  (1947).  This  method  is  used 
to  evaluate  the  back-scattering  cross  sections  for  water  spheres  at  X  =  16.230  cm  in  the  sire  range  from 

0.6  <  2to/X  <  6. 

In  Section  4  a  method  for  measuring  the  back -scattering  from  water  spheres  is  described  briefly.  Com¬ 
parison  is  made  of  the  experimental  and  the  theoretical  results. 

Another  physical  problem  of  interest  is  the  reflection  from  a  melting  snowflake  or  ice  particle.  As  a 
first  approximation,  one  may  consider,  as  the  mathematical  model  corresponding  to  this  phenomenon,  the 
reflection  from  two  concentric  spheres  with  different  complex  dielectric  factors.  Using  the  same  method 
as  for  a  single  sphere,  it  ia  possible  to  determine  a  rigorous  formal  solution  for  this  scattering.  This  ia  done 
in  Section  5. 


2-  CROSS  SECTIONS  FOR  SCATTERING  AND  ABSORPTION 
2.1.  TOTAL  CROSS  SECTIONS 

The  concept  of  scattering  and  absorption  parameters  is  found  in  many  branches  of  physics.  The  most 
commonly  used  parameters  are  those  related  to  the  total  power  scattered  and  the  total  power  absorbed  by 


an  object  in  the  field  of  an  electromagnetic  wave.  They  may  be  defined  as  follows:  Lei 

IT',  —  the  power  density  incident  upon  the  object, 

P.  »  the  total  scattered  power, 

P,  —  the  total  absorbed  power, 

P,  *■  the  total  power  removed  from  the  incident  wave 

-  (P.  +  P.)- 


Then  the  total  scattering  cross  section  is 


the  total  absorption  cross  section  b 


and  the  total  attenuation  cross  section  b 


p. 

—  p 

(2.1) 

Wi 

F\ 

(2-2) 

wt 

(2.3) 

22.  THE  BACK-SCATTERING  CROSS  SECTION 

The  back-scattering  cross  section  7  is  a  lumped  measure  of  the  ability  of  a  scattering  obstacle  to  re¬ 
radiate  energy  in  the  direction  of  the  source.  Under  the  assumption  of  a  plane  electromagnetic  wave  inci¬ 
dent  upon  the  scatterer,  this  parameter  can  he  determined  from  a  knowledge  of  the  far-aune  nattered  field. 
It  is  commonly  defined  as  follows: 

<2'4> 

where 

Wi  ~  the  power  density  in  the  plane  wave  incident  upon  the  scatterer, 

W,  »  the  power  density  of  the  far  sone  scattered  field  in  the  direction  of  the  source  and 
r  —  the  distance  from  the  scatterer  at  which  the  evaluation  of  If  .  is  made. 


It  should  be  noted  that  for  an  isotropic  scatterer,  the  back-scattering  cross  section  o  is  equal  to  the 
total  scattering  cross  section  Q,  as  defined  by  Eq.  (2.1).  In  fact,  r  may  be  defined  in  general  as  the  total 
scattering  cross  section  of  a  fictitious  isotropic  scatterer  which  scatters  energy  in  all  directum*  with  intensity 
equal  to  that  scattered  directly  back  toward  the  source  by  the  actual  scattering  object. 


Equation  (2.4)  may  be  written  in  alternative  equivalent  forms  by  writing  H  ,  and  If  ,  more  explicitly. 
Thus  IF,  b  given  formally  by 


(2.:>) 


where  4A  b  an  element  of  area  located  a  distance  r  from  the  scatterer  in  the  direction  of  the  source.  In 
addition,  dA  —  rVU,  where  dO  b  an  element  of  solid  angle  measured  from  the  scatterer.  I  sing  this 
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where  dP,/dO  is  evaluated  in  the  direction  of  the  source. 

Since  the  radiation  incident  upon  the  scattering  object  is  assumed  to  be  a  plane  wave, 

E.’ 


r. 


2f. 


(2.8) 


where  £•  is  the  magnitude  of  the  electric  field  in  the  plane  wave  and  f#  is  the  characteristic  impedance  of 
space.  Using  Eq.  (2.8),  Eq.  (2.7)  becomes 

9  ~  £,*  rftt  '  (2,9) 

**  *®T  P0*®’  ***  die  Car  cone  from  the  sea  Merer,  the  scattered  field  is  also  essentially  a  plane  wave.  There¬ 
fore,  in  terms  of  the  far-cone  scattered  field  E„ 


W. 


|E.1‘ 

2f. 


Substituting  Eqs.  (2.B)  and  (2.10)  into  Eq.  (2.4)  yields 


a  -  4rr» 


E. 


(2.10) 


(2.11) 


Finally,  if  the  far-sooe  scattered  field  is  referred  back  to  unit  range  —  i.e.,  let  -  rE.  —  then  Eq.  (2.11) 


r  -  4t 

All  of  these  equivalent  forms  are  found  in  the  literature. 


E., 


(2.12) 


2.3.  RESTRICTIONS  ON  THE  CROSS  SECTIONS 

In  the  preceding  sections,  several  cross  sections  have  been  introduced  as  convenient  parameters  for 
indicating  the  scattering  and  absorbing  properties  of  material  objects.  So  long  as  this  is  all  that  is  implied, 
tu>  diffrcsslti?*  arise.  However,  they  are  often  associated  with  the  arbitrary  picture  that  the  Poynting  vector 
is  a  unique  measure  of  the  power  flow  through  space.  The  pitfalls  of  this  incorrect  picture  are  pointed  out 
by  King  (1945).  It  should  be  noted  that  although  all  the  cross  sections  defined  have  the  dimensions  of  an 
•**■»  nooc  these  "areas”  is  in  general  equal  to  the  physical  area  of  the  object.  The  value  of  any  one  of 
the  cross  sections  depends  in  general  on  the  wavelength,  polarisation  and  angle  of  incidence  of  the  impinging 
radiation,  as  well  as  on  the  dimensions,  shape  and  material  of  the  object. 


3.  THE  THEORY  OF  THE  SCATTERING  OF  A  PLANE  WAVE  BY  A  SPHERE 
3.1.  INTRODUCTION 

The  general  problem  of  a  plane  wave  scattered  by  a  spherical  particle  was  first  studied  in  detail  by 
Mie  (1908).  He  treated  the  problem  rigorously  according  to  classical  electrodynamics.  Briefly,  his  method 
is  as  follows:  The  Maxwell  field  equations  are  written  in  spherical  coordinates.  These  equations  are  simpli¬ 
fied  by  resolving  the  electromagnetic  field  into  transverse  electric  (TE)  and  transverse  magnetic  (TM) 
waves,  and  particular  elementary  spherical -wave  solutions  of  TE  and  TM  type  are  found.  The  complete 


i.» 


solution  for  the  scattered  field  is  then  written  formally  as  the  sum  of  all  the  elementary  solutions  with  un¬ 
known  amplitude  coefficients.  The  incident  plane  wave  is  also  expressed  in  terms  of  a  linear  combination 
of  the  elementary  spherical -wave  solutions.  Finally,  the  initially  undetermined  coefficients  are  evaluated 
by  apply  ing  the  boundary  conditions  at  the  surface  of  the  sphere.  This  completes  the  formal  solution. 

The  complete  solution  to  this  problem  is  given  concisely  by  Stratton  (1941).  Following  Hansen  (1935), 
he  develops  a  set  of  orthogonal  spherical  vector  wave  functions,  each  of  which  is  a  solution  to  the  vector 
wave  equation.  The  incident  plane  wave  is  expanded  in  terms  of  these  functions,  and  the  scattered  field  is 
expressed  formally  as  a  similar  expansion  with  unknown  amplitude  coefficients.  As  in  the  Mie  solution, 
the  unknown  coefficients  are  determined  by  apply  ing  the  boundary  conditions  at  the  surface  of  the  sphere. 
The  method  of  Stratton  is  also  given  by  L.  Goldstein  (1945),  who  shows  the  relation  between  the  coefficients 
of  Stratton  and  those  of  Mie. 

In  the  following  sections  the  method  of  Stratton  is  reviewed  following  essentially  the  notation  of  Gold¬ 
stein.  The  difficulties  of  obtaining  numerical  answers  for  the  case  of  a  sphere  with  a  complex  dielectric 
constant  are  discussed,  and  it  is  shown  how  logarithmic  derivative  functions  and  the  recurrence  formulas 
due  to  Infeld  (1947)  may  be  used  to  simplify  the  computations.  Finally,  this  method  is  used  to  obtain 
quantitative  results  for  one  special  case  involving  water  spheres. 

3.2.  FORMULATION  OF  THE  PROBLEM 

Consider  a  sphere  of  radius  a,  complex  dielectric  factor  (  and  permeability  pe  which  is  isolated  in  free 
space.*  The  interior  of  the  sphere  is  called  region  I;  the  surrounding  space  is  called  region  2.  The  center 
of  the  sphere  is  chosen  as  the  origin  of  a  rectangular  coordinate  system.  A  plane  electromagnetic  wave, 
with  the  electric  field  parallel  to  the  x  axis,  is  propagated  along  the  x  axis  and  strikes  the  sphere.  Periodic 
time  dependence  of  the  form  is  assumed  4  but  in  general  will  not  be  written  explicitly.  MKS  units  are 
used  throughout.  Quantities  with  bold  letters  represent  space  vectors;  accented  quantities  in  roman  letters 
—  e.g.»  it  J.  It,  r,  ^  —  represent  unit  vectors. 

The  incident  plane  wave  may  he  expressed  as 


B.  -  i£,  - 

(3.1) 

where 

(3.2) 

is  the  wave  number  of  space  and 

It  m  w(pe*s)*  “  Y 

(3.3) 

is  the  characteristic  velocity  of  space. 

*  -  (**)-*  - 1 
k 

(3.4) 

*  This  is  a  special  case  of  the  lore  general  problem  treated  by  Mie  and  Stratton  where  the  eurrouisdmg  medium  may  aim  have 
complex  electrical  properties.  However,  the  method  ot  solution  is  identical. 

♦  Stratton  assumes  time  dependence.  Therefore,  in  compering  equations  written  hers  with  those  of  Stratton  replace  +j 
by  — i. 
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i'm.  1.  S|ihi*-ir«l  ennrdmala  lytloa. 


The  orthogonal  spherical  vector  wave  functions  are  obtained  as  follows:  The  ecalar  wave  equation  in 
epherical  coordinates  it  solved  by  separation  of  variables  to  yield 


P."  («>•  8)  (hr).  (3.5) 

Here,  P„m  (cos  8)  is  the  associated  Legendre  poly  nomial  of  the  first  kind,  nth  degree  and  mth  order;  %J*(k r) 
is  the  spherical  Bessel  function  defined  by 


-  (j^y/a+lflhr) 

(36) 

(3-7) 

Jm+i(kr)  is  the  Bessel  function  of  the  first  kind  and  half-integer  order;  H.+,w(*r)  is  the  llankel  function  of 
the  second  kind  and  half-integer  order.  The  subscripts  e  and  o  refer  to  the  even  and  odd  #  dependence. 
The  spherical  vector  wave  functions  are  defined  as 


1‘..  ~  (3-8) 

■J-  “  v  X 


(3.9) 
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<  -  Jvxm'  •  (3.10) 

Carrying  out  the  indicated  operation*  yield* 

?•"(*»•)**"** 

•■"dr 

+  - a.*  (4r)^P.- (co.  •)“*«♦  I 
r  88  »n 

p.-  (co.f)  (3.11) 

r  tin  8  co*  ' 

m'  -  x.*(kr)  P.m  (co*  #)""»♦  I 

T  ain »  co* 

-  fc" (*r)  £  P."  (co.  #)  *?*  md  $  (3.12) 

of  ■# 

*  [tit22]  *■* w  p-'  <«»•>£«♦» 

+£5^”<‘',lsp-'<“*,C"*t 

?dr.5["-‘'’(fe)l  (313> 


The  incident  plane  ware  ia  now  expanded  in  term*  of  the  spherical  vector  wave  function*.  The  rela¬ 
tion  between  the  rectangular  and  spherical  coordinate  system*  i*  readily  *een  from  Fig.  1.  Thus, 

V*  -  e-Jkrtmt  (sin  8  cos  +  cos  8  cos  d 1  -  sind$)  (3.14) 

ye-*  -  e-*  "•*  (an  8  sin  d?  +  cos#  sin  d*  +  co*d$).  (3.15) 

Since 1e-**  and  je~^"  are  solenoids]  vectors,  they  may  be  expanded  in  term*  of  and  n,  alone.  They 

are  finite  at  r  —  0;  therefore  their  expansion*  require  Bessel  function*  of  the  first  kind.  Also,  their  d  de¬ 
pendence  requires  that  m  «  1  in  the  expansions.  Finally,  comparison  of  the  even  and  the  odd  properties 
of  Eqs.  (3.12)  and  (3.13)  with  Eqs.  (3.14)  and  (3.15)  show*  that  the  formal  expansions  are 

V*  «*  -  f  (C.‘  m,!.1  +  C„*  n.u»)  (3.16) 

«-• 

-  £(C.*a.,.»  +  CVn.,.*).  (3.17) 


The  coefficients  are  evaluated  in  the  usual  way  by  applying  the  orthogonality  condition*.  For  example,  to 
evaluate  C.1,  Eq.  (3.16)  is  multiplied  scalarly  by  sin  f  m.1..1  and  the  resulting  equation  is  integrated  over 
8  and  d.  This  gives 


m,!..1  e~  *in  «  df  df 

-  /  /  E(C.‘m.„'  + 

J%  J % 


ifdfdd. 


(3.18) 


IK 


l<i«|  the  relations 


it  iit  neen  that 


Similarly, 


£  m.u.'  n,,.'  rin#d»d*  -  0 

(3.19) 

f  m.,.' •  in.,..'  rin  «<#</*«  2r|n (n  +  1  )p I*.'0 (ir)J: 

2n  +  1 

for  n'  i*  n 

for  n‘  *  n 

(3.20) 

/  /  tm.i.'e-’*'^*  rin  -  2r(—  j)M  n(n  +  l)|i. 

J  0  J  U 

(3.21) 

(3.22) 

(3.23 

(■l—l) 

(3.23) 

idcnt  plane  wave  expansions  are 

(3.26) 

B.  -  -  -  £(-/)*  f— tt;1 

»W.-|  Ln(n+l»J 

(3-27) 

Hk-  induced  necondar)  field  in  constructed  of  two  |i«rl«.  one  appltinp  inridr  I  Ik-  nplierr  and  llie  other 
«|>|>lyin£  at  all  |x>inlH  outride.  Thene  parts  are  written  an  ev|*annionn  rimilar  to  lltonc  for  the  incident  plane 
nave  hot  will*  unkmmn  amplitude  rorffiricnln.  Tlie  part  ap|*i* inp  outride  l lie  nplierr  in  referred  to  an  the 
scattered  field  and  in  indicated  by  a  niibncripl  *.  Since  lliin  field  muni  lie  repiilar  at  infinite  anil  muni  satisfy 
tlie  radiation  condition.  Ilennel  funclirinn  of  tlie  second  kind  are  required.  Tlie  Mattered  field  in  written 


!.  -  «■£<-»■[;£  +  +AW) 

1  =  ~  ^  £(-/)*  [-r"  :  1 

»».-■  Ln(»+l)J 


Tlie  fielil  inriile  tin-  n|diere  in  referml  to  an  tlie  transmitted  field  and  it  in  indicated  by  a  niilwcript  I.  Since 
lliin  field  in  finite  air  »  0.  Kennel  fum-lionn  of  the  fir*!  kind  are  m|iiirrd.  Ainu,  tlie  free  n|*a«e  count  ail  In  are 
rr|ilareri  In  tlie  com|ilr\  conntanln  of  tlie  nphere.  Tliun,  the  comjdev  wave  nnmlirr 


*i  *  «0 -i)'  ”  «  —  i 


rr|da<en  L-,  and  tlie  coan|dev  velocity 


r,  =  (iwf)  ! 
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replace*  t%.  In  Eq.  (3.29),  «,  is  the  real  effective  absolute  dielectric  constant  and  <r.  is  the  real  effective 
conductivity.  The  transmitted  field  is  now  written 

E,  -  EiE(-y)'  [  2"  +  1.1  +  yVn.,.')  (3.32) 

B,  -  -  -  £(-/)■  [  --"-rl-l  (4.W  -  >«.' n.,.1).  (3.33) 

P)  «-i  L "(n  +  1)J 


Equations  (3.28)  and  (3.29)  represent  a  formal  solution  for  the  scattered  field.  All  that  is  needed  to 
complete  the  formal  solution  is  the  evaluation  of  the  amplitude  coefficients  a.‘  and  b.'.  This  is  done  by  apply  - 
inf  the  boundary  conditions  at  the  surface  of  the  sphere.  The  boundary  conditions  at  r  «  a  are  * 


t  X  (E.  +  E.)  -  ?  X  Bi 
?X(B,  +  B,)-^X  B,. 


The  d  component  of  Eq.  (3.34)  yields 

P, 


dP.a)  (cos  0)  [  (os.'”  («)!'  +  (a)l'  b.'lNa*^  {Na)Y 


+  j 


[ 


fra 


4iO 


0. 


The  0  component  of  Eq.  (3.34)  yields 

| +  o. *»(«*)  - a.'*.(,)(JV«) J 


+j 


.  P.1  (cos  *)  I  {a*,<0  (a))'  +  6„,[a*,<,)  (a)]'  6,/|iVas«(l)  (.Vo)!' 


(co«<)  I 

sin  0  \ 


ka 


*P 


0. 


The  d  component  of  Eq.  (3.35)  yields 


P„*  (cos 0)  (*.«>(<*)  +  A.'z,'11  (a)  bm‘x.w0\a)  ) 

*  r,  ) 


ml 


dfV  (cost)  f  |m.(l> (a)]'  +  0.101. w(a)r  _  o.'l \W" (/Vo))' 

r»*i« 


vjte 

The  0  component  of  Eq.  (3.35)  yields 

dP.1  (cos  0 )  |  s.(l>(o)  +6.’s.t!:(a)  _  t.‘s.<l>(Aa)  [ 
40  I  tv  -  til 


0. 


+  j 


.  P.1  (cos  0) 


|os,<,)(a)|'  +  «.*[a*.,’,(«)l'  n^W"(^Vu)|' 


vjca 


ijna 


0. 


(3.34) 

(3.35) 


(3.36) 


(3.37) 


(3.38) 


(3.39) 


Equations  (3.36),  (3.37),  (3.38)  and  (2.39)  are  all  satisfied  if  the  quantities  in  the  brace*  |  |  are  identically 
equal  to  aero.  This  yields 

«.'*.<0((Va)  -  «.**.“*(«)  -  *."’(*>  <»•»> 


-  Equation  (MS)  takes  this  timplii  fans,  without  the  magnetic  constant*,  Mitre  lash  rrfiia  I  end  region  2  err  wessd  to  beta 
the  same  permeability. 
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a.1iV«s.‘“(iV«)r  -  e.1«.*(«)r  -  [«.<*>(«)]' 

(3.41) 

Nb.,s.*)(No)  -  *.■*.»(«)  -  *.«><«) 

(3.42) 

MNTa».®(iV«)r  -  Nb.1cm. *(«)]'  -  JVl«s.w  («)]'. 

(3.43) 

These  rimuha 

aeons  equations  are  solved  for  the  nnudficients.  Thus, 

(NaM«*®(«)r  -  s.w  (*,)[JW°  (Na)Y 
s.01  (Ne)[«x.®  [a)Y  ~  *.»  («)(AT«.W  ( Na))’ 

(3.44) 

.  «.w  (•)(*•«.<»  (Na)Y  -  N*m.<°  (JVo)[as.«»  («))' 

s.w(n)lAr«s.w(JVn)|'  -  N*s»®  (Na)[ax.®  (a))* 

(3.45) 

Here  the  prime*  at  the  equate  bracket*  indicate  difleremiatiot)  with  respect  to  the  argument  of  the  Bessel 
hactws  inside  the  brackets,  a  ■  km  —  2rs/X;  Na  —  io;  N  ■  ky Jk  ■  ({/*»)*  is  the  complex  index  of  re* 
fraction  for  the  sphere. 


3.3.  THE  FAR-ZONE  FIELD  AND  TOTAL  POWER  RELATIONS 

To  determine  the  scattered  and  absorbed  power,  it  is  necessary  to  know  the  Car-eone  scattered  field  — 
the  field  farr>a  This  is  found  by  introducing  the  asymptotic  value  for  the  spherical  Hankd  function 

(hr)  *  ^  •-**'-  (*+ (3.46) 

kr 


into  the  expansion  formulas  for  the  scattered  field,  Eqs.  (3-28)  and  (329),  and  neglecting  terms  of  higher 
order  than  (1/r).  When  this  is  done,  it  is  seen  that  the  far-eooe  scattered  field  in  component  form  is 


Car  ■  Bit  ■  0 

(3.47) 

(3.4«) 

(3.49) 

The  total  field  at  any  point  external  to  the  sphere  is  obtained  ky  vector  addition  of  the  incident  and 

scattered  fields. 

Thus, 

»  -  +  B.;  B  -  B,  +  B.. 

(320) 

The  complex  Poynting  vector  corresponding  to  this  total  field  is  defined  by 

8- jIXB*. 

2 

(321) 

The  asterisk  ind 

Beats*  Um  ooaplex  conjugate,  and  *  is  the  Magnetic  constant  of  apace  (*  ■  1/^). 

Since 

the  far-aone  field  is  transverse,  the  Poyating  vector  is  radiaL  Substituting  Eq.  (3.50)  into  Eq.  (3.51)  and 
carrying  out  the  vector  product  yields 
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S  -  fS,  -  ?[j|  (E..BS  -  -  E^B, 


+  ~2  (E,#S  +  E*BJ  -  E+B.,*  -  •  (3.52) 

The  couples  Pointing  vector  corresponding  to  the  incident  field  is 

S.  -  9S„  «  ?  ~  (EitB+*  -  E*Bn*),  (3.53) 

end  that  corresponding  to  the  scattered  field  is 

S.  -  iS„  «  (E,tBS  -  £„£!.,*).  (3..VI) 

The  total  scattered  power  is  obtained  by  integrating  S„  over  the  surface  of  a  large  sphere  of  radius  r 
and  taking  the  real  part.  Thus, 


P.-\  R« 


n: 


(E.,B+*  -  )H  sin  #d»W*. 


The  total  power  absorbed  by  the  spltere  is  obtained  by  integrating  S,  over  the  surface  of  a  large  sphere  and 
taking  the  negative  of  the  real  part.  Thus, 


-  Jit* 


ff 


S jr*  sin  1 4 


The  total  power  removed  from  the  incident  wave  is  equal  to  the  sum  of  P.  and  P,.  Denote  this  by  Pi. 
Since  the  integral  of  Sir  over  a  closed  surface  is  aero,  P,  can  be  expressed  as  an  integral  of  the  third  term  in 
Eq.  (3.52).  Thus, 

P,  «  P.  +  P,  -  -  ^  r*  sin  I S  d*.  (3.57) 

The  integral  for  P,  is  evaluated  as  follows:  From  Eqs.  (3.48)  and  (3.49) 

E.&S--(— )  £  £  ([—-+. IT2**  1  1 

«%\kr/  ,_i  lLn("  +  1)J  b*»(«*  +  1)J 

r  .  - p-p-  .  L  ^  ^ ip-* *p-'  .  .*  p.‘  rf#*.1 

•  n-’o.**  .  -  +  *.•*.**  —  —  +  «.*b.**  —  —— 

L  sus1*  4  4  sin  t  4 

jp  i  p  n  ) 

+ *-■■*--*  ~ir  ifij  * }  <“•' 

and 

_  ns  f  ir ± n  rj*±jLi  r  ,  - gy gv 

E-  ®*  ■ « w.?1„?1ll«(« + dJ  L-(-+»J  L--  *  * 

Mil2  I  OT  Nil  $ 

P  1  >~|  ) 

+  fc"lB"'*5fi_dfJain’4r  (VW) 
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Substitute  Eqs.  (3.58)  and  (3,59)  into  Eq.  (3.55)  and  recall  that 


Then, 


r 


CO a**d* 


sin*  d  M  -  »• 


P.  - 


tmR1 

2a** 


-JT 

.?1^TT5]Lot] 

..  /  p,'  dp,'  (  dp.»  p.'  \ 

"  \ain  I  4  dl  sin  1/ 

.  .  „/dP.'P..'  P.«  dPjYl  .  1 

Ur  ST#  +  dT#  -dT;J  - ']  *■ 


+ 


Using  the  following  identities, 

for  ail  n  and  m 

J %  \su>  9  49  49  sin  9/ 


r/dp,'dpm'pm'pj\  . 

J*\  49  49  sin1  #  /  * 


949 


0  fern's 

2[«<n  +  1)P 


(2«  +  1) 


for  n  ■»  as. 


one  gets 


By  a 


process, 


P.  -  ^ZT  £  (2«  +  1)  (o.*n.-  +  6.V). 

•%**  ..i 


Pi  -  -  Re  £  (2n  +  !)(«.*  +  6.‘). 


(3.60) 

(3.61) 

(3.62) 

(3.63) 

(3.64) 


P,  and  P<  may  be  eapraaaed  in  terms  of  the  cross  sections  defined  in  Section  2.  Thus,  using  Eq.  (2.1), 
the  total  scattering  cross  section  is 

Q.  -  5  £  (2n  +  l)(n.*n.-  +  6.6.-).  (3.65) 

r  ..t 

and,  using  Eq.  (23),  the  total  attenuation  cross  section  (including  both  absorption  and  scattering)  ia 

Q,  -  -  ^  Re  £  (2n  +  l)(o.*  +  6.*).  (3.66) 

*r  .-I 


3.4.  THE  BACK-SCATTERINC  CROSS  SECTION 


It  is  now  possible  to  determine  the  back-scattering  cross  section  defined  in  Section  2.2.  Frans  Eq.  (2.9) 


Frans  Eq.  (3.55) 


8rf.  dP, 

£.*  da  t, 


E+  Bo*)- 


(3.67) 

(*M) 
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Unag  Eqe.  (3.48)  aad  (3.49),  aad  noting  that 

_  ^dfVtom#)  _  (  i)«*if,>("  +  1}T 
•-»»  #-»r  a  L  2  j 

Hi* MM  that 

~fj  -  ^  t  L  (-l)'+“(2*.  +  1)(2m  +  !)(•.■  -  *.*)(•.**  -  *.-)• 

•U  p«v  ••  1  **l 

HmnCdn, 

#  -  n  £  £  (-i)’+-(2a  +  l)(2m  +  1)(«.*  -  *.*)(«.-  -  k.**). 

*  «-l 


(3.W) 

(3.70) 

(3.71) 


Aa  ahacftatire  method  of  deriving  the  equation  for  the  back  scattering  croaa  aection  ia  to  uw  the  equiva- 
efimtion  at  a  inroiriag  only  the  for-uone  Mattered  field.  In  view  of  the  orientation  of  the  coordinate 
ay  Mem  need  (aee  Fig.  1).  Eq.  (Ml)  ia 


E-C--.n> 


(3.72) 


9  ■  4wr*\ 


^Cv.) 


Uarag  Eqk  (3.48)  aad  (3.49)  together  with  Eq.  (3.69)  yielda 

'EiLriEti-  (t)|.?,(-TJ)(-  ’)'(-• -v)i 


♦-»/* 


-  (s)  l£(b +,)(- ')'(-• 


(3.73) 


(3.74) 


(3.75) 


Equation  (3.75)  ie  idMtical  with  Eq.  (3.71). 


3.5.  THE  SCATTERING  AMPLITUDE  COEFFICIENTS 

From  the  [—**■;  eactlono  it  ia  eeea  that  the  evaluation  of  the  far-aone  aeatteted  field  and  the  virioua 
croea  eeetione  redneea  to  the  determination  of  the  ecatteriag  amplitude  coefficient!, •„*  aad  I/,  and  the  auroma- 
tion  of  the  appropriate  aeriea.  The  aeattered  field  may  be  regarded  aa  due  to  the  forced  oactHatfem  of  mag¬ 
netic  aad  electric  mahipolct.  If  thia  point  of  new  ia  taken,  it  ia  eaay  to  aee  that  the  coefficienta  «,*  are 
aaaociated  with  eeciBatieaa  of  atagaetic  type,  aad  the  coefficient*  b»*  are  aaaociated  with  oariBarion*  of  electric 
type.  It  ia  recalled  (Eq.  (3.12))  that  the  rector  ware  function  n'^  hae  no  radial  component.  Therefore,  if 

the  coefficient*  b.*  are  al  aero  aad  the  toainau  a.*  are  not  aero,  the  R  field  ia  purely  tmaarene  aad  the 
B  field  hae  a  radial  component.  (Son  Eqe.  (3.28)  aad  (3£9).)  Siaulariy,  if  the  coefficfoata  a.*  are  all  aero 
aad  the  cotfficitnl*  k»*  are  aot.  the  E.  field  hm  a  radial  component  aad  the  B,  field  ie  purely  traaereroe.  The 
magnetic  aad  electric  typee  of  ooriRation  are  aemetfaam  called  traaereree  electric  (TE)  aad  traaerene 
magnetic  (TM)  typee,  respectively. 

Aa  Mated  above,  the  coefficient*  a.*  aad  i,'  amy  be  eeaoeiated  with  the  forced  oeciletione  of  nugnetic 
aad  electric  mahipolm.  Whenever  the  forced  frequency  approaches  one  of  the  natural  frequeadee  of  vibra¬ 
tion,  a  condition  of  reeMaare  occur*.  The  natural  freqncarice  of  ribra  tiau  am  determined  by  set  ting  the 
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denoniiv.  .f>-„  equal  to  zero  in  Eq».  (3.44)  and  (3.45).  However,  the  natural  frequencies  are  complex  while 
the  impressed  frequency  is  always  real.  Therefore,  no  difficulties  arise  at  resonance. 

The  problem  of  the  numerical  evaluation  of  the  scattering  amplitude  coefficients  will  now  be  considered. 
This  problem  is  solved  formallt  by  Kqs.  (3.44)  and  (3.45).  However,  they  are  quite  formidable,  and  atten¬ 
tion  is  usually  directed  toward  the  special  cases  where  certain  simplifications  are  possible.  These  will  he 
reviewed,  following  Stratton  (1941),  before  the  more  difficult  general  problem  is  taken  up. 

For  a  metal  sphere,  |/Var|  >  1  and  the  asymptotic  expressions  for  the  Bessel  functions  may  be  used. 
They  are 

x.01  (iVa)  *  ^  cos  [nc  -  — ^  1)Tj ,  (3.76) 


lNea.{'\Na)]'  *  -  sin  ^Na  -  —  j 


Substituting  these  expressions  into  Eqs.  (3.44)  snd  (3.45)  yields 


U(,>(n)| 

u««r 

([«.U)(a)rl 

1  [«.w  (a)]'} 


The  final  equations  are  exact  for  a  perfectly  conducting  sphere. 

For  a  large  sphere,  a  >  1,  the  asymptotic  expressions  may  be  used  for  all  the  Bessel  functions  involved 
in  Kqs.  (3.44)  and  (3.45).  The  result  is 

.  .  ,  ...  ^  r«in  A  —  /V  cos  A  Ian  A#"] 

L  l+.w..„«  J  •  <SJ»> 


.  .  /  ...  *  r**n  A  —  /V  cos  A  Ian  Af 

‘  ("y)^L  I  +  //V  tan  Af . 

,  .  <+1  j.  ("cos  A  tan  Af  —  yV  sin  A"| 

*'  *1 - tan  Af  —  /TV - J' 


cos  A  un  Af  — 


(»  +  0« 


Af  =  Va  - 


("  +  1)t 


Tliese  equations  are  included  because  they  illustrate  the  general  oscillatory  behavior  of  the  scattering  ampli¬ 
tude  coefficients. 

For  a  small  sphere,  s<  1,  it  is  possible  to  get  a  good  approximation  for  the  scattering  amplitude  coeffi¬ 
cients  by  substituting  the  series  expansions  of  the  spherical  Bessel  and  Hankel  functions  into  Eqs.  (3.44) 
and  (3.45)  and  keeping  only  a  few  terms  of  the  resulting  se-irs.  This  has  been  done  by  L.  Goldstein  (1945), 
who  corrects  the  results  of  Stratton  (1941).  The  series  expansions  of  the  spherical  functions  are 


_<i>  f  \  „  o’ a"  T  f  ( _ II  I 

{»/  *  r< -p"(»  +  w)!  i 

*.  (■)  -  “  .^tLm!(2n  +  2m  +  l)!j< 

/  y.  f*  (2w  —  2«)ll 

+  21«'+,.t.Lm!(n-m)!j^ 


Substituting  these  expansions  into  Eqs.  (3.44)  and  (3.45)  and  keeping  only  a  few  terms, 

•MEriJ-kss-JI--] 


•[l  +  o5 


(n/V>  +  n  +  1) 

(2n  +  1 )[ (2r»  -  1)IV»  -  n  -  1| 
(2n  +  3)(2n  -  l)(niV»  +  n  +  1) 


(3.85) 


_  ~  )’  .j(2»  +  l)(n  +  l)(y_^l_)l  +  . .  1  (3.86) 

7  (2n  +  1)!  J  l  (nAf*  +  n  +  1)  j  J 


If  a  is  so  small  that  powers  of  a  higher  than  the  sixth  may  be  neglected,  only  three  coefficients  are  needed: 


«i*  * 


zl 

45 


((V*  -  1)b* 


(3.87) 


(3.88) 

(3.89) 


This  is  an  important  case  since  it  applies  to  the  raindrop  problem  for  many  of  the  radars  in  use.  Equations 
(3.87),  (3.88)  and  (3.89)  are  the  same  ones  used  by  Ryde  (1941,  1944)  if  account  is  taken  of  the  difference  in 
definition  between  the  coefficients  defined  by  Stratton  (used  here)  and  those  of  Mie  (uaed  by  Ryde).  The 
relation  between  the  two  sets  of  definitions  is  given  by  Goldstein  (1945): 

p»M**  -  (-  l)"/(2n  +  1)«.*  (3.90) 

c.**  -  (-  l)*+1/(2n  +  l)b.‘.  (3.91) 

It  should  be  noted  that  if  a  is  so  small  that  powers  of  a  higher  than  the  third  may  be  neglected,  only  the 
electric  dipole  mode  is  needed.  This  is  the  case  of  Rayleigh  scattering. 

It  is  seen  that  even  with  the  simplifications  afforded  by  the  approximations  involved  in  the  special  cases, 
the  numerical  evaluation  of  the  scattering  amplitude  coefficients  is  quite  involved.  An  exact  evaluation 
would  be  even  more  so.  Goldstein  (1945)  remarks,  "An  exact  computation  of  these  coefficients  is  out  of  the 
question  on  account  of  the  lack  of  tables  of  Bessel  and  Hankel  functions  of  complex  argument  needed  here. 
The  special  cases  discussed,  however,  do  not  cover  all  the  problems  of  interest,  so  that  it  is  necessary  to  con¬ 
sider  the  problem  of  getting  reliable  answers  in  the  general  case. 

Although  the  evaluation  of  the  scattering  amplitude  coefficients  is  difficult  in  the  general  case,  it  can  be 
accomplished.  In  their  study  of  the  effects  of  meteorological  elements  on  microwave  radiation,  Ryde  and 
Ryde  (1945)  set  up  a  schedule  for  computing  the  coefficients  directly  by  calculating  the  values  of  the  spherical 
Bessel  and  Hankel  functions  for  the  arguments  needed  and  substituting  them  into  the  Mie  equivalent  of 
Eqs.  (3.44)  and  (3.45).  Since  the  work  by  Ryde  and  Ryde,  Lowan  (1949)  has  computed  the  coefficients  for 
various  values  of  a  and  six  values  of  JV  corresponding  to  water  in  the  microwave  region.  There  are  still  no 
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extensive  tablet  available,  however,  for  the  ipherical  functions  of  complex  arguments,  and  there  is  no  general 
rwiImmI  for  getting  exact  values  of  the  coefficients  except  the  "brute  force”  method.  In  Section  3.6,  such  a 
meilKxl  Hill  lie  given. 


3.6.  EVALl  ATION  OF  THE  SCATTERING  AMPLITUDE  COEFFICIENTS  USING  LOGARITH- 
MIC  DERIVATIVE  FUNCTIONS 


The  general  solution  for  the  scattering  amplitude  coefficients  is  given  b>  Eqs.  (3.44)  and  (3.45). 
These  equations  max  be  rearranged  as  follows: 


(3.92) 


(3.93) 


It  should  be  noted  that  in  this  form  the  equations  involve  the  logarithmic  derivatives  of  {ox,(l>  (a)], 
[/Vax„0)(./Va)|  and  [a».l,,(a)j.  Logarithmic  derivatives  of  this  type  were  used  bv  Infeld  (1947)  in  his  study 
of  the  spherical  antenna  with  a  gap.  and  hi  Smith  (1918)  and  Tai  (1949)  in  the  study  of  bk-onical  antennas. 
If  the  notation  of  Smith  and  Tai  is  used,  the  Bessel  density  functions  are  ilefined 


S.<*)  =  .(*)],  (3.94) 

K.(x)  -  |r1//.+1(!,(r)].  (3.95) 


Following  Tai  (1949),  tlie  logarithmic  derivatives  of  tliese  Bessel  density  functions  are  defined  as  new  func¬ 
tions,  c.(x)  and  *,(*).  Thus, 


*.(*) 


(3.%) 


.  JTM  I 

P'[  K.(x)  “  [**.<*(*)]  ' 


(3.97) 


In  terms  of  these  functions,  the  scattering  amplitude  coefficients  are 

.  _  _  [*»(«)  -  N*m{Na)- 1 

“  Ul,,(a)Jlp.(a)  -N<r.(JVa)J 


b.‘ 


_  r^)l  fr.(Na)  -  Ne.(«)1 

L*.w(o)J’U(No)  -7VP.(o)J 


(3.98) 

(3.99) 


This  farm  appears  to  be  simpler  since  there  are  fewer  functions  and  no  derivatives  of  functions  involved 
explicitly.  To  make  the  apparent  simplification  real,  however,  it  is  necessary  to  demonstrate  that  values  for 
the  logarithmic  derivative  functions  can  be  easily  obtained.  This  will  now  be  done. 

Since  the  following  derivations  hold  for  all  the  Besael  density  functions,  involving  either  the  Bessel, 
Neumann  or  llankel  functions,  it  is  convenient  to  let 


•>" 


[/.<*)] 

Z,(x)  —  •  \.(x )  1  =  the  general  evlindrical  function, 

«.<*) 

(.3.100) 

D.  (x )  ■  x*Z,+|(x)  “  the  general  Hessel  densilv  function  and 

(3.101) 

/  v 

!«(*)  ■  q  ^  ^  m  the  general  logarithmic  derivative  function. 

(3.102) 

Differentiate  Eq.  (3.101)  with  respect  to  x: 

&.(*)  -  *‘Z'.+!(x)  +  |x-‘Z.+,(x). 

(3.103) 

Substitute  the  recurrence  formula  * 

7'n+.(x)  -  Z._.(x)  -  [^]z.+!(x) 

(3.101) 

into  Kqi  (3.103)  to  eliminate  tlie  derivative  of  Z: 

iy.(x)  =  x:  -  ”z.+!(x)J- 

• 

(3.105) 

Divide  bv  D,(x): 

,  M ,  iy-M.  ,  r«.-L«i  «. 

D«(x)  L^»4:U)-I  ■*' 

(.3.106) 

This  is  the  standard  equation  used  to  evaluate  tlie  logarithmic  derivative  functions  wlicn  x  if-  real. 


The  difficulty  comes  when  x  is  complex.  Kq  nation  (3.106)  cannot  lie  used  then  because  tables  of 


Z»+|(x)  are  not  available  for  x  complex.  In  this  case,  use  is  mailr  of  tlie  recurrence 
(1947),  which  will  now  be  derived.  From  Kq.  (3.101), 

formula  doc  to  Infeld, 

Z„+j(.v)  -  x~ilh(x). 

Differentiate: 

(3.107) 

z',+.(*)  />.(*)!- 

Differentiate  again: 

(3.106) 

Z"«+j(x)  -  x-s[D".(x)  -  ~  IF.Cc)  +  (JA  lh(x)\. 

The  liesael  equation  of  half-integer  order  is 

(3.109) 

*-*£".+.(*)  +  xZ',+i(x)  +  [x-  -  (n  +  J)s|Z„4i(jr)  -  0. 

Substitute  Eqs.  (3.107),  (3.108)  ami  (3.100)  into  Kq.  (3.110)  and  collrii  terms. 

(3.110) 

n"Jx)  =[^t  “  -  ']/»•(*)• 

(3.1 II) 

This  is  tlie  same  as  Kg.  (D.l)  of  Infeld  if  n  is  replaced  bv  p  —  }.  It  is  now  necessarv  to  pet  an  expres¬ 
sion  for  D'.(x).  Substitute  Kqs,  (3.10T)  and  (3.106)  into  the  recurrence  furi.ula* 


Sn,  tviuapb,  l  itna.  A  Trtmtiwe  am  the  7V»;i  sf  Btmtti  Faam'sai,  pp.  45  and  66. 
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Z.+1(*)  -  f11— *]z„_,(*)  -  Z'._,(*) 

(3.112) 

and  collect  terms: 

/>.(*)  -*Dm.dx)  -P'._,(*)- 

(3.113) 

Differentiate: 

p'«(*>  -  -  p'.-tw  -  2  ®-iw  -  »".-.(*)• 

X  Xr 

(3.114) 

Substitute  Eq.  (3.111)  into  Eq.  (3.114)  and  collect  terms: 

P'.W  -  -  O'..  ,(*>  -  5  D.-ito  +  (3.115) 

X  Xr 

Now,  divide  Eq.  (3.115)  by  Eq.  (3.113).  Then, 


But, 


Therefore, 


«.(*) 


PM*) 

P.(*) 


5zy._l(*)-^D._,(*)  +  D._,(*) 


5z)._t(*)  -  2y..,(*) 


*«-i(*)- 


«.(*) 


=  !._,(,)  -=j+l 
j|  ~  «.-!(*) 


(3.116) 


(3.117) 


This  is  the  same  as  Eq.  (D.6)  of  Infold  if  n  is  replaced  by  p  —  J.  It  it  sometimes  easier  to  use  Eq.  (3.117) 
if  both  numerator  and  denominator  are  cleared  of  fractions.  Thus, 


x*  +  nx4,_i(s)  —  n* 
nx  -  x*4._,(x) 


(3.118) 


Equation  (3.118)  gives  the  recurrence  relation  whereby  the  logarithmic  derivative  function  of  any 
order  can  be  found  if  the  function  of  the  next  lower  order  is  known.  It  will  now  be  shown  that  the  lowest 
order  can  be  computed  directly  .  Since  c.(x)  it  the  only  function  with  complex  argument  in  Eqs.  (3.98) 
and  (3.99),  only  this  function  will  be  considered: 


$•(*)  -  l**/t(*)]  -  «n  *. 

(3.119) 

s'.(*)  -  k‘/»(*)r  -  (-Vccx. 

v/ 

(3.120) 

Therefore, 

W  S»(*) 

(3.121) 

Now,  if  x  -  c  —  jd,  then. 

/  \  «  «n  2c  +  /  ainh  2d 

**w  Skui'^.u 

(3.122) 

Higher -order  terms  may  now  be  found  using  Eq.  (3.118). 

It  i*  seen  that  by  using  Eqs.  (3.98)  end  (3.99)  together  with  Eqs.  (3.106),  (3.118)  and  (3.122),  the 
scattering  amplitude  coefficients  of  any  order  can  be  computed  exactlx,  even  in  regions  where  the  spherical 
Bessel  functions  of  complex  arguments  are  not  tabulated.  In  the  next  section  this  method  is  applied  to  the 
evaluation  of  the  hack -scattering  cross  section  for  one  particular  rase  involving  water  spheres  with  sises 
comparable  to  the  wavelength.  This  case  has  also  been  studied  expert menlalU  and  comparison  between 
the  theoretical  and  experimental  results  is  made  in  Section  4. 


3.7.  NUMERICAL  COMPUTATIONS 

The  method  given  in  Section  3.6  waa  used  to  compute  the  scattering  amplitmle  coefficients  for  water 
spheres  at  X  ■»  16.230  cm  and  0.6  <  a  <  6.  These  coefficients  were  then  used  to  compute  the  back -scattering 
cross  section,  using  Eq.  (3.75).  For  plotting  purposes,  it  is  convenient  to  normalise  the  luck-mattering 
crass  section  with  respect  to  the  geometrical  cross  section.  If  this  is  done,  Eq.  (3.75)  becomes 

(3.123) 

This  is  the  equation  actually  used  in  the  computations.  Since  the  values  of  o«*  and  6.*  decrease  rapidly 
for  n  >  «*,  it  is  necessary  to  carry  the  summation  only  to  n  ^  2a. 

Since  the  computed  values  of  the  back -scattering  cross  section  depend  on  the  index  of  refraction,  it  is 
important  to  have  reliable  values  for  this  quantity.  It  is  recalled  that  with  m  *>  pj,  the  complex  index  of 
refraction  is  related  to  the  complex  dielectric  factor  by 

n  -  0‘  -  <d*  -  Vr  -  yro‘.  (3.124) 

where  {r  is  the  relative  dielectric  factor,  and  and  —  f",  are  the  real  and  imaginary  parts  of  (r.  Ryde 
(1941,  1944)  and  Goldstein  (1945)  determined  by  using  the  Debye  (1929)  formula. 


9 

WO* 


2 

a* 


£  (2n  +  1)  (  -  1)*  («.*  -  *.*) 

»-i 


+ 


~  **  , 
l+/(VX) 


(3.125) 


where  is  the  relative  optical  dielectric  factor;  $,»  i*  the  relative  static  dielectric  factor;  and  X«  is  the 
transition  wavelength.  They  determined  X*  from  Eq.  (3.125)  by  using  Collie's  (1944)  measured  values  of 
and  ("r  for  X  ■  1.26  cm.  The  validity  of  using  the  Debye  equations  for  the  frequency  dependence  of  £, 
at  constant  temperature  has  been  verified  by  Saxton  (1945)  and  Collie  et  al.  (1948).  Slight  differences 
have  been  found  in  the  experimental  values  of  X*,  but  in  general  the  agreement  is  good.  This  same  method 
has  accordingly  been  used  in  the  present  research.  Hie  result  is  f,  *  81  —  /7.8  at  X  ■  16.230  and  temper¬ 
atures  near  20s  C. 


Some  of  the  computed  values  for  the  logarithmic  derivative  functions  and  scattering  amplitude  coeffi¬ 
cients  are  plotted  in  Figs.  2-19. 

Figures  2  and  3  show  *„(a)  versus  a  for  n  —  1,2, 3, ... ,  12.  The  values  of  v. (a)  were  computed  using 
Eq.  (3.106)  with  Z„±t(a)  —  /„A ,(<*).  It  is  seen  that  r,(a)  is  a  smooth,  non-oscillating  function.  It  ap¬ 
proaches  +  <*>  as  (n  +  i)/a  when  a  approaches  aero.  It  has  poles  at  the  aeros  of  /m+t(a).  The  poles  of 
r.(a),  however,  do  not  cause  the  scattering  amplitude  coefficients  to  blow  up  since  it  is  always  the  product 
/„+l(a) -#.(a)  which  appear*  in  Eqs.  (3.98)  and  (3.99),  and  this  product  has  no  poles. 
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Km;.  IS.  TW  mllrrini  amplituiie  nnSriwO,  h  ’  ami  fc<*.  Km:.  19.  The  waUerimf  amplilwle  rarlKiMiK  h’.  k1  ami  k'. 

Figure*  4  and  5  allow  the  real  and  imaginary  part*  of  p.  (a).  The  value*  of  these  function*  were  com¬ 
puted  using  Eq.  (3.106)  with  Zm±i(a)  ■  W,±|W(a).  It  i*  seen  that  they  are  also  non -oscillatory.  The  real 
part  of  p.(a)  approaches  -  *  ai  —  m/m  when  a  approaches  aero.  A*  a  increases,  the  real  part  increase* 
rapidly  from  —  »  and  approach**  aero.  The  imaginary  part  starts  at  aero  and  approaches  —  1  as  a  increases. 

Figures  6  to  13  dhow  the  real  and  imaginary  parts  of  »(Afo)  for  n  —  1, 2, . . . ,  8.  It  is  seen  that  these 
are  rapidly  oscillating,  highly  damped  functions.  The  real  part  approaches  aero  and  the  imaginary  part 
approaches  I  a*  a  increases. 

Figure*  14  to  19  show  the  scattering  amplitude  coefficients  a.*  and  &.*,  for  n  —  1,  2, ....  7.  The  general 
oscillatory  behavior  of  the  coefficients  is  seen,  as  well  as  the  rapid  decrease  in  amplitude  for  n  >  a. 

The  theoretical  normalised  back  ■scattering  cross  section  curve  is  shown  in  Fig.  20,  where  it  is  compared 
with  the  one  obtained  experimentally. 

4.  EXPERIMENTAL  DETERMINATION  OF  THE  BACK-SCATTERING  FROM 
WATER  SPHERES  AND  COMPARISON  WITH  THEORY 

4.1.  INTRODUCTION 

Since  the  introduction  of  microwave  techniques  lo  the  field  of  meteorology,  much  work  has  been  done 
on  the  specific  applications  to  the  detection  of  reflections  from  rain  and  other  hydromeleora.  Unfortunately, 
most  of  this  work  has  been  qualitative  in  nature.  Some  quantitative  measurements  on  rain  have  been 


nude  by  H.  Goldstein  (1945),  Marshall  et  aL  (1947),  Lanplle  and  Gunn  (1948)  and  Hooper  and  Kippax 
(1950).  rwnfcrini;  the  number  of  variable*  involved,  them  measurements  show  good  agreement  with 
theory.  It  is  difficult,  however,  to  isolate  the  individual  effects  when  there  are  many  variables,  and  it 
seems  desirable  to  have  some  laboratory  measurements  on  single  spheres.  Recently,  this  need  has  been  at 
least  partially  filled  by  the  work  of  Klotebaugh  and  Duckett  (1949)  on  plastic  sphere*  simulating  raindrops 
and  the  work  of  Aden  (1950)  on  actual  water  spheres.  In  Sections  4.2  to  4.4,  the  techniques  used  by  Aden 
will  be  described  briefly,  and  the  results  obtained  will  be  compared  with  those  determined  from  theory. 

A2.  THE  EXPERIMENTAL  PROBLEM 

The  iiipisimniilal  measurement  of  the  reflection  of  electromagnetic  waves  from  individual  water  spheres 
involves  two  important  difficulties:  (1)  Since  the  sphere  is  a  very  low  gain  ref.ee  lor.  it  is  hard  to  obtain 
reliable  answers  by  ordinary  poising  techniques  and  (2)  it  is  not  easy  to  maintain  a  water  sphere  while 
measurements  are  being  made.  For  these  resMons,  no  direct  laboratory  measurements  had  been  under¬ 
taken  successfully  prior  to  the  work  mentioned  above.  Even  now,  these  difficulties  almost  preclude  the 
taking  of  measurements  in  the  region  of  greatest  interests  —  be.,  the  spectral  region  where  water  spheres  of 
actual  raindrop  aise  are  comparable  to  the  wavelength.  However,  it  is  possible  to  perform  an  experiment  at 
slightly  longer  wavelengths  and  to  increase  the  sise  of  the  water  spheres  to  get  into  the  critical  sire  region. 
Ibis  is  the  experiment  that  actually  was  performed.  It  is  believed  that  this  technique  is  useful  ainoe  com¬ 
parison  between  experiment  and  theory  at  one  frequency  should  give  a  good  indication  of  the  correlation 
to  be  expected  at  other  frequencies. 

A3.  THE  EXPERIMENTAL  TECHNIQUES 

The  experimental  method  need  to  measure  the  back -scat  taring  cross  section  was  the  standing-wave 
method  of  D.  King  (1948).  This  method  utilises  an  image  screen  technique,  and  the  standing  waves  set 
up  on  the  screen  by  the  interaction  of  the  incident  and  reradiated  waves  are  measured  along  the  radial  line 
between  the  scatterer  and  the  source.  This  method  offers  the  advantages  of  a  system  having  absolute 
calibration,  relatively  simple  equipment  and  measurements  at  low  power.  The  main  disadvantages  are 
that  it  requires  an  image  screen  which  is  large,  uniform  and  rigid,  and  that  obtaining  results  is  very  time 
consuming.  The  derivation  of  the  formulas  needed  and  a  diecuarion  of  the  approximations,  as  well  as  a  de¬ 
tailed  description  of  the  equipment  used,  are  given  by  Aden  (1950). 

One  of  the  worst  «h«»srlw  to  ■«»*fc«ng  measurements  on  water  spheres  is  the  inability  to  maintain  such 
spheres  while  measurements  are  being  taken.  This  obstacle  was  overcome  by  using  as  a  container  thin 
hnmii|iluii  it  il  shell  forma  of  Styrofoam*  mounted  on  aluminum  disks.  Since  Styrofoam  has  dielectric 
properties  extremely  dose  to  thorn  of  air,  it  had  a  negligible  effect  on  the  measurements.  Thus,  when  a 
form  was  filled  with  water  and  the  disk  was  inserted  into  its  proper  place  in  the  image  screen,  the  effect 
wee  that  of  a  hemisphere  of  water  exposed  over  a  large  image  screen.  By  image  theory,  the  measurements 
taken  were  the  same  as  for  i  conplcte  iplwt  it  free  ipacc- 

4.4.  COMPARISON  OF  THE  THEORETICAL  AND  EXPERIMENTAL  RESULTS 

Measurements  of  the  back -scattering  crate  section  were  made  on  30  water  spheres  in  the  elec  tries  I  sise 
region  0.74  <  a  <  5.90.  In  Fig.  20,  the  results  of  these  measurements  are  plotted,  together  with  those 
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determined  from  theory.  The  theoretical  curve  was  obtained  by  using  the  method  of  logarithmic  deriva¬ 
tive  functions  given  in  Section  3.  It  is  seen  that  there  is  very  good  agreement  between  the  experimental 
and  theoretical  results.  For  spheres  very  much  smaller  in  sue  than  those  considered  here,  the  results  should 
be  given  to  s  good  approximation  by  the  Rayleigh  law.  For  spheres  of  larger  siae  than  those  considered 
here,  the  results  should  approach  the  large -else  approximation  of  about  0.64. 


5.  THE  THEORY  OF  THE  SCATTERING  OF  A  PLANE  WAVE  BY  A  SPHERE 
WITH  A  CONCENTRIC  SPHERICAL  SHELL 

5.1.  INTRODUCTION 

In  the  preceding  sections,  the  scattering  of  electromagnetic  radiation  by  a  single  sphere  has  been  con¬ 
sidered,  and  its  application  to  the  more  general  problem  of  microwave  reflection  from  rain  has  been  indicated. 
Another  physical  problem  of  interest  is  the  scattering  of  microwave  radiation  by  melting  snow  and  ice 
particles.  For  purposes  of  analysis,  this  physical  problem  may,  to  a  first  approximation,  be  replaced  by  the 
mathematical  problem  of  finding  the  reflection  from  a  sphere  with  a  concentric  spherical  shell  of  different 
dielectric  factor.  This  Utter  problem  can  be  solved  rigorously. 

The  method  of  solution  is  a  direct  extension  of  that  used  for  the  problem  of  the  scattering  from  a  single 
sphere.  The  incident  pUne  wave  is  expanded  in  terms  of  the  orthogonal  spherical  vector  wave  functions  of 
Stratton  (1941).  The  induced  secondary  fields  in  the  various  regions  are  written  as  similar  expansions 
with  unknown  amplitude  coefficients.  As  in  the  previous  case,  the  unknown  coefficients  are  determined 
from  the  boundary  conditions.  In  this  case,  however,  the  matching  must  be  done  simultaneously  over 
two  boundary  surfaces  instead  of  one. 

5.2.  FORMULATION  AND  SOLUTION  OF  THE  PROBLEM 

Consider  a  sphere  of  radius  a,  complex  dielectric  factor  £i  and  permeability  p*  which  is  surrounded  bv 
a  spherical  shell  of  inner  radius  o  and  outer  radius  b  with  a  complex  dielectric  factor  fi  and  permeability  /j*. 
The  resulting  configuration  is  assumed  to  be  isolated  in  free  space.  The  interior  of  the  sphere,  the  interior 
of  the  shell  and  the  surrounding  space  are  called  regions  1, 2  and  3,  respectively.  As  in  Section  2,  the  center 
of  the  sphere  is  chosen  as  the  origin  of  a  rectangular  coordinate  system;  the  incident  electromagnetic  wave  is 
propagated  along  the  s  axis,  and  the  electric  vector  is  linearly  polarised  parallel  to  the  x  direction.  This  is 
illustrated  in  Fig.  21. 

With  the  conditions  stated  above,  the  expressions  for  the  incident  plane  wave  are  exactly  the  same  as 
in  Section  3: 

-  \E.  -  tE*-*  -  E,  £  f  2"  +  )  1  (m.,.1  +  /*„.')  (5.1) 
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6  *%  a-1  T  t;  j 

The  induced  secondary  field  must  now  be  constructed  in  three  parts,  one  applying  in  each  of  the  three 
regions  defined  above.  These  parts  ate  written  as  expansions  to  those  for  the  incident  wave  but 
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Fie.  21 .  Plaar  »*ve  imidrnt  upon  a  sphere  with  a  concentric  inimical  shell. 

with  unknown  amplitude  coefficient*.  The  part  applying  outside  the  shell  is  again  referred  to  as  the  scat¬ 
tered  field  and  ia  indicated  by  a  subscript  i.  The  form  of  this  expansion  is  the  same  as  in  Eqs.  (3.28)  and 


(3.29)  although  the  values  of  a*  and  &'  are  different  here: 

E  (  —  j)*  F  ~7  T  i\l  (•■‘'■bu*  n,!.’)  (5.3) 

■-i  Lnvn  + 

"•  -  -  f  .?,< (5*> 

By  analogy  with  the  previous  problem,  the  field  inside  the  sphere  is  called  the  transmitted  field,  and  it  is 
indicated  by  a  subscript  t.  Again,  the  formal  expansions  are  the  same  as  before  (Eqs.  (3.32)  and  (3.33)) 
with  the  understanding  that  a,1  and  I.1  have  different  values: 

£,-£.£  (-  jy  [w2("  +  («.'mel„'  -1- jbm%lm')  (5.5) 

*  e.  •  *  r  2*  tin 

B*  -  -  r  E  (  - »"  -T-I  n  -  ya.'n.,.1).  (5.6) 

**  «-i  Ln(»  +  1)J 

Here,  t>i  —  (mfi)-*  is  the  complex  characteristic  velocity  of  the  sphere. 

Inside  the  spherical  shell,  the  terms  involving  Bessel  functions  of  both  first  and  second  kinds  must  be 
retained.  The  formal  expansions  may  be  written 

*.-££(-/)’  +  d.’m,,.1  -I- /(bA.i.*  +  B.X  ..*))  (5.7) 

■*-“,£(- /)' [;^  + p] +  B.’m,,.*  - /(s.^.'  +  /I.W)]  (541) 
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when  •»’.  A,\  k,*,  £„’  are  unknown  amplitude  coefficient*,  and  •*  ■  (m^)  *  i*  the  couple*  characteristic 
Telocity  of  the  ebeiL 

Equation*  (5.3)  through  (5.8)  represent  a  formal  eohition  far  the  induced  secondary  field.  AD  that  i* 
needed  to  complete  the  formal  eoiutiofi  i*  the  erahutioa  of  the  eight  amplitude  coefficient*.  Thi*  i*  done 
by  applying  the  boundary  condition*  at  the  two  surface*  of  dielectric  discontinuity.  The  boundary  condi¬ 
tion*  at  r  -  a  are 


tXE,|  -*X£.| 

r-i  ml 

(5.9) 

9  X  B«|  -?Xfit| 

(5.10) 

and  at  r  ■  b  are 

9X(E,  +  E.)|  -  9  X  Ei| 

r-J  r-t 

(5.11) 

9x(B,  +  B.)|  -fxB,|  . 

r-»  r«» 

(5.12) 

Theee  lead  to  two  *et*  of  simultaneous  equation*  each  involving  four  unknown*  a*  follow*: 

[JW >(JV,a)r 

Ni 
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-*.®(»)fi.*  +  0 
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and 
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-  *.*(»)  +  o 
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-  [**.*<-)]'  *.*  +  0 
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Here,  a  »  2m/X ,  r  —  2 *4/X  ;  JV|  —  VfJ/i  and  iVi  -  Vfa/n  are  the  complex  indice*  of  refraction  for 
the  ephere  and  shell,  respectively.  A*  before,  the  prime*  at  the  aqnare  bracket*  indicate  differentiation  with 
respect  to  the  argument  of  the  B easel  function  inside  the  bracket*.  These  two  set*  of  aimultaneou*  equa¬ 
tion#  may  he  solved  for  the  eight  amplitude  coefficient*. 

To  evaluate  the  back-acattering  cross  section,  it  i*  necessary  only  to  solve  for  the  scattering  amplitude 
coefficient*, «,'  and  i.'.  If  this  is  done,  and  if  the  derivative*  of  the  Bessel  density  functions  are  eliminated 
by  introducing  the  logarithmic  derivative  functions,  the  result  is: 

*.‘°(r)  ["Ft  +  ••(rJCi’l 
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The  equation  for  the  bark-neat  terinp  *  row  net-lion  in  tlw  name  an  before: 
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Kquatkm  (5.27).  loprllier  with  Kqn.  (5.21)  tlimuph  (5.26).  pi  nth  the  formal  solution  to  tlie  problem  of 
back  M-atterinp  from  a  njilterr  with  a  t-onrrntrit-  n|ilierit-al  nliell.*  \s  a  rltrrL  it  should  lie  noted  that 
thin  nolution  should  leilwe  to  that  for  a  sinple  n|>liere  ife  ■  i  anti  V,  —  V.  ■  V.  Sinee  the  equation  for 
the  Itat-k-niatln-inp  crow  seel  ion  in  the  name  in  hotli  t-anen.  thin  rrquiren  the  neatlerinp  amplitude  coeffirients 
t«i  lie  tlie  name.  Kxamination  of  Kqs.  (5.21)  through  (5.26)  reveals  that  ibex  ilo  rrdtire  to  Kqs.  (3.68) 
anti  (3.66)  under  the  conditions  ntaleil. 

Mllioitph  tlie  a|i|iliealitin  of  thin  |irobleiii  to  |utrlit'iilar  plnniral  problems  han  tieeu  inilieatrtl.  the  detailn 
of  tliene  apidiealionn  will  ihm  In-  given  here.  Tliev  an-,  however,  lieing  pivrn  attention  liv  other  memlters 
of  tlie  (•rtqilivnit'n  lleneart'lt  Division. 


Sill*  ■-  lit**  *■* ini |il< *1  miii  til  ilii*  iii.iii«i>i  rt|il  4ihI  kitkit  Imic  Invlnl  tlir  bhmt  grthtal  riNr  itlhir  nil  llimr  rrtgiaim 
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